Rheology of cholesteric blue phases 
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Blue phases of cholesteric liquid crystals offer a spectacular example of naturally occurring discli- 
nation line networks. Here we numerically solve the hydrodynamic equations of motion to investigate 
the response of three types of blue phases to an imposed Poiseuille flow. We show that shear forces 
bend and twist and can unzip the disclination lines. Under gentle forcing the network opposes 
the flow and the apparent viscosity is significantly higher than that of an isotropic liquid. With 
increased forcing we find strong shear thinning corresponding to the disruption of the defect net- 
work. As the viscosity starts to drop, the imposed flow sets the network into motion. Disclinations 
break-up and re-form with their neighbours in the flow direction. This gives rise to oscillations in 
the time-dependent measurement of the average stress. 
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In many cholesteric liquid crystals the transition be- 
tween the regular cholesteric phase and the isotropic 
phase occurs through one or more intermediate phases 
known as "blue phases" (BPs) [1-3]. Understanding the 
structure of the blue phases was a theoretical and exper- 
imental tour de force clearly summarised in [3] . 

Blue phases occur because the locally favoured state 
is a double twist (i.e. a rotation of the director field 
in two directions, as opposed to the single twist of the 
cholesteric) which is globally incompatible with the re- 
quirement of continuity [3] . This results in the formation 
of networks of defects that separate local regions of dou- 
ble twist with axes aligned in different directions. 

The networks can be intricate and lead to structures 
that are periodic at large length scales. For example 
the BPI and BPII phases are 3-dimensional cubic phases 
where the orientational order can be periodic on the scale 
of the wavelength of visible light. Due to this, they have 
potential applications in fast light modulators [4], tun- 
able photonic crystals [5] and lasers [6] . The main factor 
limiting the technological exploitation of BPs is that they 
are generally only stable over a temperature range ~ 1 K. 
However recently it has been shown that the addition of 
polymers can stabilise BPI over more than 60 K opening 
perspectives in BP technology. 

Though the static structure of BPs is well accepted es- 
sentially nothing is known of their rheology. Given the 
intricacy of the director configurations one might expect 
that elastic and viscous responses to an external stress 
combine to give highly non-Newtonian behaviour under 
flow. Cholesterics show a strong dependence of viscos- 
ity on stress and, in particular, a striking permeation 
mode where the backflow acts to minimise director dis- 
tortions [8]. In blue phases there is also a fully three- 
dimensional network of defects whose elastic properties 
will increase the richness of their rheology. 

Therefore this paper aims to investigate the rheology of 
cholesteric BPs by numerically solving the hydrodynamic 



equations of motion for BPI, BPII and, for comparison, a 
simplified metastable network which we will term double 
twist (DT). We consider a Poiseuille flow field. Under 
gentle forcing we find that the viscosity is consistently 
higher than that of an equivalent isotropic phase but, as 
the forcing is increased, there is strong shear thinning 
corresponding to the progressive rupture of the network. 
As the disclination network starts to flow the periodic- 
ity of the blue phase gives rise to an oscillating state 
which would leave a measurable signature in the time- 
dependent stress curve. Elastic networks are of interest 
in many contexts and we compare our results to very re- 
cent studies of the rheological properties of disclinations 
in colloidal suspensions in liquid crystals [9]. 

The equilibrium properties of BPs are described by 
a Landau-de Gennes free energy density, written in 
terms of a tensor order parameter Qap [1,3]. This com- 
prises a bulk term 
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and a distortion term [3] 
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where K is the elastic constant and p = 27r/g'o is the 
pitch of the cholesteric. Aq is a constant, 7 controls the 
magnitude of order and e^cs is the Levi-Civita tensor. 

The time evolution of the system is governed by the 
equation of motion 
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where F is a collective rotational diffusion constant and 
Dt is the material derivative for rod-like molecules [10]. 
The molecular field H^p is explicitly given by 
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The fluid velocity, iT, obeys the usual continuity equation 
and the Navier- Stokes equation 
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The stress tensor Ila/3 is 
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where p is the fluid density, 77 is an isotropic viscosity and 
^ is related to the aspect ratio of the molecules. Pq is a 
constant in the simulations reported here. Note that the 
order parameter field affects the dynamics of the flow field 
through the stress tensor. This is the backfiow coupling. 
Details of the equation of motion can be found in [8] . To 
solve these equations we use a three dimensional lattice 
Boltzmann algorithm [11,12]. 

In the first column of Fig. 1 we show the equilibrium 
configuration of the three different BPs we shall consider. 
The cylinders are disclination lines between regions of 
double twist which are repeated periodically to build up 
the blue phase structure. Within the tubes the order 
parameter attains a deep local minimum, whose value is 
slightly flow-dependent. Figs. IB and IC are the accepted 
structures of BPI and BPII. The configuration shown in 
lA, is the DT structure. It is the simplest disclination line 
network that can be constructed with double twist, i.e. a 
periodic array of doubly twisted cylinders with axes along 
z. DT is met ast able and has not been observed. However 
its simplicity makes it a useful tool to help understand 
the rheology of the other blue phases. Local regions of 
the DT configuration are likely to be found during a tem- 
perature quench while a third observed blue phase BPIII 
is thought to comprise a gas of doubly twisted cylinders 
like the ones building up the DT network [1,3]. 

Each structure was obtained by relaxing to its min- 
imum by solving Eq. 3 numerically. Periodic boundary 
conditions were used and initial configurations were set 
according to the approximate solutions in [3]. For more 
details and an equilibrium phase diagram see [13]. 

To study the rheology of BPs, we sandwiched them 
between two plates a distance L = p/2 apart along the 
2;- axis. (Qualitatively similar results are obtained with 
larger L.) Poiseuille flow along y was imposed by a pres- 
sure gradient AP together with no-slip conditions for 
the velocity at the plates. We quantify the strength of 



the forcing via the dimensionless parameter F = , 
where c is the sound velocity. At the boundaries we as- 
sume that the director profile corresponding to the BP 
structures does not change under flow, that is the discli- 
nations are fixed at the boundaries. Physically, this is 
likely since small surface irregularities should pin the 
equilibrium director profile. A similar situation occurs 
in permeation [8], where fixed boundary conditions give 
results closer to the experiments. 

Fig. 2 summarises the dependence of the apparent 
viscosity, r^app, on the strength of the forcing and com- 
pares the blue phase behaviour to that of a corresponding 
Newtonian isotropic fluid, r^app is found by computing 
the flux of u through a plane perpendicular to y and by 
comparing it with that of a Newtonian fluid in the same 
geometry and under the same body force. For small forc- 
ing, F ^ 10~^, r^app increases by a factor of ~ 4 times 
over that of the isotropic fluid. All three phases reach 
a stationary state in which the networks of line defects 
is deformed by the flow, as shown in column II of Fig. 
1. The deformations consist of first bending, and then 
twisting, the disclination line network without changing 
its topology. In BP I in particular we observe that the 
disclinations twist either separately or in pairs of different 
handedness, as predicted in [1,2]. 

Backfiow is significant and acts to decrease the bend- 
ing and twisting of the network lines. This can be clearly 
seen by comparing the results of the full simulations in 
column II in Fig. 1, to those in column III, obtained 
in calculations in which backfiow was neglected and the 
order parameter was simply advected by the flow. 

The primary flow is symmetric to a good approxima- 
tion for BPI and BPII but displays an asymmetry for DT. 
This is due to the different orientation of the disclination 
loop with respect to the flow near the top and bottom 
plates, which arises as the two disclinations twist. In the 
case of BPI and BPII the disclinations are more numer- 
ous, so that any asymmetry is washed out by averaging. 
In the steady state we find non-negligible secondary flows 
along both the x and z directions. 

For small forcing the viscosities of BPI and BPII do 
not vary greatly. The DT phase, however, shows shear 
thicknening. This is because of a spectacular change in 
the defect structure in which each line of 27r disclinations 
(defects of topological strength 1) characterising the DT 
structure opens up to a disclination ring (defect of topo- 
logical strength 1/2). The ring then bends and twists on 
itself as the flow increases. A similar transition, driven 
by applying an external magnetic field, has been studied 
in a system of colloids dispersed in a nematic [14]. 

We now consider what happens as F is increased. Fig. 
2 gives evidence of significant shear thinning for all three 
structures considered. The disclination network is de- 
stroyed by the flow and the structure becomes that of a 
standard nematic subject to a Poiseuille flow [1,2]. 

The high and low viscosity behaviours are separated 
by a crossover region in which the viscosity decreases 
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rapidly. This corresponds to an intriguing network dy- 
namics, shown, for the doubly twisted structure, in Fig. 
3. The disclination lines in a given cell break up and 
part of the defect structure moves with the flow into the 
neighbouring cell. It then re-forms the network by at- 
taching to the portions of the defect left behind. The 
process repeats and hence there is a continuous process 
of destruction and rebuilding of the textures. The dy- 
namics can be monitored through the time dependence 
of the free energy and the stress tensor of the system, 
both of which oscillate (Fig. 4), at least up to the 1 ms 
time scale reached in the simulations. 

Early experiments have reported high viscosities in 
BPs [15]. It is also interesting to compare our results 
with those obtained very recently in Refs. [9] for a discli- 
nation line network obtained by dispersing spherical col- 
loidal particles in a cholesteric liquid crystal. These au- 
thors measured the apparent viscosity versus shear rate 
and found a shear thinning behaviour that closely resem- 
bles that shown in Fig. 2. Finally we comment on the 
rheology of cholest erics. Here there is an even higher vis- 
cosity at low shear rate, now because the director field 
(rather than a disclination network) does not want to 
distort. In both systems there is, with fixed boundary 
conditions, strong shear thinning at higher shear rates 
[8] . In cholesterics when fiow is perpendicular to the he- 
lix a nematic-like viscosity is observed [16]. 

In conclusion, we have reported the first numerical in- 
vestigation of the rheology of blue phases in cholesterics. 
We applied Poiseuille flow starting from three topolog- 
ically different equilibrium defect networks: two corre- 
spond to the accepted structures of BPI and BPII, and 
the other is a periodic array of doubly twisted cylinders. 
Under small forcing, the network opposes the flow giv- 
ing rise to a signiflcant increase in the apparent viscosity 
77app. Shear forces cause the disclinations to bend and 
twist. In the case of the doubly twisted structure we 
observe a flow- induced unzipping of the array of lines 
into a network of loops which causes shear thickening at 
small forcing. As the forcing is increased, we find clear 
evidence of shear thinning. Similar behaviour is seen in 
rheological experiments on colloidal suspensions in liquid 
crystals. In the crossover region, as r^app drops, we also 
predict a novel oscillatory regime which results from the 
interplay of the periodicity of the blue phases with the 
imposed fiow. Here the network continuously breaks and 
re-forms as portions of the disclinations in the centre of 
the channel move to neighbouring cells and relink with 
the parts of the network left behind by the fiow. 

We have considered a director field pinned on the 
boundaries at z = 0,p/2. Experience with permeative 
fiows [8] suggests that if free boundary conditions are 
used there should be no increase in r^app, although the 
rest of the phenomenology should be preserved. 

Our algorithm should be applicable to study the rheol- 
ogy of other 3-dimensional disclination networks in liquid 
crystals or in liquid crystalline gels, such as those sta- 
bilised by colloids [9], emulsions or polymers. An inter- 



esting related problem is the switching in BP devices [7] . 
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FIG. 1. Equilibrium and dynamic steady states for 
cholesteric BPs. Each row corresponds to a different BP: A. 
DT; B. BPI; C. BPII. Each column corresponds to a different 
flow configuration: L no flow; 11. Poiseuille flow with back- 
flow; III. Poiseuille flow without backflow. Configurations are 
drawn in terms of their disclination networks, each tube sur- 
rounds a disclination and marks where the order parameter 
drops to 60 % of its maximum value (the threshold determines 
the width of the tubes, e.g. if it is chosen as 50 % the width 
becomes ~ 10% smaller). The director field is doubly twisted 
between disclinations. Parameters were K = 2.5 10~^^ N, 
L — p/2 — 0.5 /xm, 7 = 3, while F ^ 3.3 10~^ for columns II 
and III; while the Leslie viscosity ratio a2>loLi ^ 0.08, so the 
BP was flow aligning. BPs are often described by a pair of 
dimensionless number, their chirality k and reduced temper- 
ature r (defined in [3,4]). In our case = 0.59 and r = 0.35. 
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FIG. 2. Apparent viscosities of BP I, BPII, DT, and of a 
reference isotropic fluid (I) , as a function oi F — . 
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FIG. 3. Time evolution of a DT structure under Poiseuille 
flow, for F — 0.023, for which the system is in the crossover 
regime, see Fig. 2. 2x2x1 structural periods are shown and 
distinguished by different colours (the actual structure is infi- 
nite along both x and y). The disclinations first bend (frame 
2) and then unzip (3). The unzipped structure near z = L/2 
is then extended in the direction of the fiow (4) and merges 
with its neighbour along y (5). There is some relaxation (6) 
and then the centres of the disclination are again stretched by 
the flow field (7,8) until the structure (9) resembles (4) and 
the process repeats. Parameters are as in Fig. 1. The times 
corresponding to each of the frame are also shown. 
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FIG. 4. Averaged stress (solid line) and free energy density 
(dashed line) corresponding to the dynamic evolution shown 
in Fig. 3. Oscillations are apparent in both curves. 
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